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Abstract 

For Belavin's elliptic quantum R-matrix, we construct an L-operator as a set of difference operators acting 
on functions on the type A weight space. According to the fundamental relation RLL — LLR, the trace of 
the L-operator gives a commuting difference operators. We show that for the above mentioned L-operator 
this approach gives Macdonald type operators with elliptic theta function coefficient, actually equivalent to 
Ruijsenaars' operators. The relationship between the difference L-operator and Krichever's Lax matrix is 
given, and an explicit formula for elliptic commuting differential operators is derived. We also study the 
invariant subspace for the system which is spanned by symmetric theta functions on the weight space. 



1 Introduction 



In [Ml], [M2], I. G. Macdonald defined a commuting system of difference operators for each root system and a 
new family of orthogonal polynomials containing two rational parameters (q, t) (in case all the roots have the 
equal length) as their similtcncous eigenfunctions. Up to now, there are at least two theories which provide 
understanding of the system with general value of (q, t). One is the work by Etingof and Kirillov, who obtained 
these operators as the image of central elements of the quantum enveloping algebra U q {sl n ) acting on "vector 
valued characters" [EK1|. The other is the work by Cherednik |C92|, who used double affine Hecke algebra, its 
representation via g-difference operators, and the center of the algebra. (See also Section 7.) 
Here we wish to suggest yet another approach for the system. 

Needless to say, the Yang-Baxter equation is one of the important backgrounds of the above two works. 
Originally, in Baxter's study of two-dimensional lattice statistical models, the Yang-Baxter equation arose as 
the condition to provide sufficiently many commuting operators. This is done by taking the trace of the so-called 
L-operators, the operators which satisfy the " RLL = LLR relation" (^). In the lattice model situation, the 
R-matrix reads as the local Boltzmann weight of the model. Given an R-matrix satisfying the Yang-Baxter 
equation, an L-o perator natu rally arises as the row-to-row transfer matrix and then its trace gives the commuting 
transfer matrix |Bax71 [TF|. 

Apart from the lattice models, recall the Lax matrix method in the completely integrable systems. Once 
the equation in problem is formulated in a Lax form, natural candidates of commuting integrals of the motion 
are the characteristic polynomials of the Lax matrix, namely the traces of the powers of the matrix. Their 
commutativity does not hold in general but one can make use of the r-matrix structure, a differential form of 
the RLL = LLR, relation, to ensure it (see |S93| and the references therein). 

Therefore the following question seems quite natural to ask: "what kind of operator arise if we start with the 
L-operator realized as difference operators for appropriate functions and take the trace?" This is our approach 
and we will show that this idea actually works quite well at least for one interesting case. 

The case we consider in this paper is for the elliptic R-matrix of Belavin [Be] and the main goal is the 
reproduction of Ruijesenaars' elliptic system. In the trigonometric limit, up to a certain simple "gauge trans- 
formation" [Re] , this R-matrix degenerates to the image of the universal R-matrix for the quantum affine 



enveloping algebra U q (A n _ 1 ) (Q,[Dr[) in the vector representation. The n = 2 case of Belavin's R-matrix 

*Short title : Ruijsenaars' operators as commuting transfer matrices. AMS subject classification: 81R50, 17B37, 05E05. 
T Partly supported by The Kawai Foundation for the Promotion of Mathematical Science, Japan Society for the Promotion of 
Science, and Grant-in-Aid for Scientific Research on Priority Areas, the Ministry of Education, Science, Sports and Culture, Japan. 
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is nothing but Baxter's eight vertex model, for which the novel quadratic algebra is defined by Sklyanin |Hj. 
We hope that the present paper will give some insight into the structure related to the Sklyanin algebra and 
Belavin's elliptic R-matrix. 

The plan of this paper is as follows. For Belavin's n-state R-matrix, we have constructed an L-operator 



whose matrix elements are a certain difference operators in the previous paper [HI | (Section 2, Theorem 1). This 
is an sl n - generalization of the L-operator corresponding to the Sklyanin's difference operators in [Bj. For this 
L-operator, we apply the fusion procedure (Section 3) and compute the traces of the fused L-operator (Section 
4). The RLL = LLR relation ensures the commutativity of the resulting operators. It turns out that the 
traces exactly give us Macdonald-type difference operators Md (d = 1, • • • , n), whose coefficients are given in 
terms of Jacobi's elliptic theta function (Theorem 2, main result). Fixing the elliptic modulus parameter r, 
Belavin's R-matrix has one parameter h, while the factorized L-operator admits another parameter c as well as 
the spectral parameter u. Consequently, the operator Md depends on these three parameters h, c, u and it turns 
out that the spectral parameter u appears only in the overall factor. The remaining two parameters just play 
the role of Macdonald's parameters q,t and the explicit correspondance is given by q = cxp TT\/—lfi, t — q~^ . 

The computation of the trace of the fused L-operators hinges upon an interesting theta function determinant 
formula (Lemma 1). It contains the parameter h. In the limit h — > 0, this formula degenerates to the well-known 



Cauchy type determinant, which probably first appeared as formula (12) of [Frob] in the literature, also known 
as the genus one case of Fay's trisecant formula fFayfl . 

In Section 5, we give explicit relationship between our elliptic commuting operators and other approaches 
for the system. First of all, consider the trigonometric limit. Then it is easy to see that our operators tend 
to Macdonald's operators. In this trigonometric case another presentation of the system is known, that is, 
the trigonometric Ruijsenaars' operators. The equivalence of Macdonald's operators and Ruijsenaars' ones is 



given via the conjugation by a certain function multiplication (due to T.Koornwinder; see [Die]). In Subsection 
5.1 (Proposition 2) we extend this equivalence to the present elliptic case, namely between our commuting 
difference operators and Ruijsenaars' relativistic elliptic Calogero-Moser system M. Secondly, we look at the 



differential limit case, i.e. the elliptic Calogero-Moser system. This system admits the Lax formalism [Kr|. Now 
the question is whether and how Krichever's Lax matrix arise from the factorized L-operator; this is answered 
in Subsection 5.2 (Proposition 3). Subsection 5.3 is devoted to the generating function for the operators {Md}- 
Since Md is defined to be the trace of the factorized L-operator in the "degree d exterior" fused representation, 
we expect that the generating function is just the characteristic polynomial of the L-operator. This is in fact 
the case: Theorem 3. The formula can be considered as an elliptic extension of the generating function in 
Macdonald's case (Q,VI(3.2)). An analogue of Jiro Sekiguchi's generating function [pel is obtained as well. 
Further, it turns out that our formula is useful to derive the commuting differential operators in an explicit way 



we will give an elliptic generalization of A.Debiard's De] formula 



As in the theory of Macdonald polynomials, one may consider the family of symmetric functions defined 
as the joint eigenfunction for {Md}d=i,--,n- This diagonalization problem is still under investigation. Here 
we would like to establish some structure theorem on a certain invariant subspace for the L-operator. Setting 
the parameter c to be a nonnegative integer I, it is the space Thf spanned by the level I A^-i affine Lie 
algebra characters, or the symmetric theta functions on the weight space of type A. This space is a higher rank 
analogue of Sklyanin's finite dimensional function space. We will state that this space can be identified with 
the symmetrically fused representation as the module of the L-operator algebra (Theorem 5) . 

Brief discussion is given in Section 7. Appendix is devoted to proving the theta function identity Lemma 1. 



2 Review of the factorized L operator 

For n > 1 let V = © fce Z/nZ Ce ' £ ^ = ek+ ") and let 9> h e GL(V)be ge k := e k exp^,he k := e k+1 . We have 
gh = hgexp^ 1 . Let h, r € C, Imr > 0. We assume Ti £ Z + Zr for convenience. 

Belavin's R-matrix R(u) = R%(u) is characterized as the unique solution of the following five conditions. 

• Rn(u) is a holomorphic End(V <X> V) -valued function in u, 

• Rn{u) = (x <g> x)Rn(u)(x ® x)^ 1 for x — g,h, 

• R n {u + 1) = (, 9 ® l)- 1 ^^)^ <g> 1) x (-1), 

• R r (u + t) = (h^^Rniu)^®!)- 1 x (-exp27ri(u + | + 
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Rn(0) = P ■ x ® y >—> y <8> x. 



One verifies that 1) there is a unique solution to the above conditions and 2) the solution satisfies the Yang- 
Baxter equation. An explicit formula for R(u) is also available. Put 



i{u,t) := ^2 cxp27ri(/iu+ — r) 

fi£m-\-l'Zi 



(i) 



and 0W(it) := 6i_£ i(u + \-, tit). The zeroes of 9^\u) are given by Z + (j + nZ)r. Then we have [RT] 



i',j' = l 



(2) 



See figure 1 for our convention to visualize the matrix R. 

[ Fig. 1 ] 

By an L-operator we mean a matrix L(u) — [L(v,)j]i t j = i t ... n of operators (noncommutative letters) that 
satisfies 



R(u - v)L(u) (g> L(v) = L(v) (g) L(u)R(u - v), 



(3) 



where R(u) := PR{u). 

For Belavin's R-matrix we shall construct such an L-operator in the following way. Let h* be the weight 
space for sl n (C). We realize h* in C n = ©»=i,... jn Cej, < e^Cj >= S^j, as the orthogonal complement to 
Si=i ■■■ n ei - We denote the orthogonal projection of by e^. For each A, \x £ h* and j = 1, • • • , n we define the 
so-called intertwining vectors ([ Bax73| for n = 2 case and generalized in JMOl ) 



9j(%- < A, efc >)/v = Tr?(r) : /x - A = ?Lefc for some fc = 1, • 
: otherwise 



, n. 



(4) 



where 



,•(«) 



M(u+ 2 )+ 2^ T 



and r/(r) := p 1 / 24 IIm=i — P m ) denotes the Dedekind eta function with p := exp 2ttit. 

[ Fig. 2 ] 

Let further define (j>{u) l t +hek ' : > to be the entry in the inverse matrix to [4>( u )a +hek j]j,k=i,---.n, namely [QF] 



£^(«C l8 * j >(«){: +RE *'i = 'v..A. 



(5) 



fe=i 



The pictorial expression for (/>, <fi are given in FigsJ^J^ respectively, and then (|^) can be written as in Fig.|]. 

[ Fig. 3 ] 



[ Fig. 4 ] 

Then generalizing a result in the paper |s| (II) , we have 
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Theorem 1 ([Hl|, [H2|) For a function f on h* , put 



(L(c|u)j/)(/*) := £ 4>(u + chr+^j^ur+^fif, + he k ) 



(6) 



fc=i 



('Fig.^j. Then for any c e C, </ie collection of difference operators L(c\u) = [L(c\u)j]i l j=i... in satisfies the desired 
relation (Qj. i.e., L(c\u) gives a 1-parameter (c) family of ' L- operators. 



[ Fig. 5 ] 



Remark. This type of "factorized" L-operator first appeared in [[K|. It was for the trigonometric two by 
two R-matrix case, namely the six vertex model corresponding to the algebra U q (sl2). Later, in connection with 
the chiral-Potts and the Kashiwara-Miwa solutions of the Yang-Baxter equation, generalization is given in this 
form in the U q (sl n ) case by [BKMS] and the elliptic case by [QF] and |H1|. 

Let us give a b rief acc ount for the proof of this theorem. In the proof, the so-called face model 

Boltzmann weight [ JMQ1 | (abbrev. face weight) plays a central role. The face weight is a collection of the 
quantities (Figure fig:ffacewt) 



W 



A u v 



[ Fig. 6 ] 

defined for A, /i, \i! ', i/£h' and u e C. Explicitly they are given by the formulas 

_ 9(u + h) 



W 



A + hei 
X u A + 2hi~i 
A + Mi 



W 



W 



A + he~i 
A u A + h(e~i + e~j) 
A + he~i 

A + he~i 
A u A + h(e~i + <Tj) 
A + he-j 



8(h) ' 

_ ej-u + Xjj) 
eiXij) 

9(u) 6(h + Xij) 
9(h) 9(X i3 ) 



(i + 3), 
(< * 3), 



where :=< A, — ej >, and 

oo 

9(u) := 9i A (u + -) = V^lp 1/S (z 1/2 ~ z- 1 ' 2 ) TJ (1 - Z p m )(l - z^p m )(l - p rn ) 



(7) 



(z 1 / 2 — exp tts/— lit, p = exp 2tt^/— It) denotes the Jacobi theta function. For the other configurations of A, /i, p! 
and v, the face weight is set to 0: 

A* 

W A u v :=0. 
The face weight and Belavin's R matrix are related by the following relation (Fig.fr]), 



£ r(u ~ vt^mti ® = E ® m^vft 



A u — v v 



(8) 
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[ Fig. 7 ] 



i.e. they are "intertwined" by the intertwining vectors. It follows that the face weight satisfies the so-called 
face Yang-Baxter equation (see (pjl))- By the duality relation (ph, (0) immediately implies 



E t(vr/ 



(u)fR(u-v)y., = y w 



(9) 



[ Fig. 8 ] 

(Fig. ||) which is in fact equivalent to (||) itself. Combining these two relations, we have 



A* 

A u — v v 



(Fig. 



i',j';tM' 

||) and the theorem follows. 



Fig. 9 ] 



3 Fusion Procedure 



3.1 Overview 

Let be the bialgebra generated by the formal letters {L{u)^\u E C,i,j = 1, • • • ,n} with respect to the 

relation (||) and let C(h*) be the ring of meromorphic functions on h*. Recall V = ®j=i,--- tn Ce^ ~ C™. Our 
L-operator (^]), giving a representation L(u) l j i— > L(c\u)j of A(R) by the definition, is an endomorphism on the 
space V® <D(h*) : 

L(c\u) G End(T/®0(h*)). 

Here the first space V = C" should be regarded as the defining comodule (vector "co" representation) for A(R) 
: the comodule structure depends on the spectral parameter u and is given by 



U u ■ V 3 e l 



n 

7 = 1 



®L{u)) €V®A(R). 



We denote the space V endowed with this comodule structure by V(B lt ). We can consider more complicated 
comodules for A(R) as well. That is, for each Young diagram Y and a parameter u we can construct an 
yl(i?)-comodule V(Y U ) whose dimension is just the same as for the GL„-module that corresponds to Y. In 
contrast with the trigonometric R-matrix case, where we have the quantized enveloping algebra UniA^^) and 
its universal R-matrix as the origin of the tensor category structure, in the elliptic R-matrix case such an 
algebraized theory does not seem to be available until now. Yet we can utilize a constructive definition of 
V(Y U ) known as the fusion technique | KRS ||C|. This is done by taking appropriate sub/quotient of the tensor 
comodule V(Q Ul ® • ■ • ®Uu k ) = V(Q U1 ) <£> • • • <£> V(U Uh ) ~ V® k . Applying this fusion technique to our case, we 
get a collection of difference operators 



L(c\Y u ) G End(VTO®0(h*)) 
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for each Y, as well as the "fused R-matrices" R(Y U ,Y^) : V(Y U ) <8 V(YJ) -> V(Y„') <8 V(Ya). The latter gives an 
isomorphism between the two A(i?)-comodules V(Y U ) ® V(Y„') and V{Yv) ® ^C^u) an d it follows that we have 
the relation 

fl(y U) rj) l(c|f u ) ® L(c|r t ;) = l(c|y„') o l( c |y») fl(y u , y„'). (10) 
3.2 Antisymmetric case 

For our purpose of constructing Macdonald type commuting difference operators, important are the fused L 
operators corresponding to the Young diagrams of vertical boxes that we are now going to describe. 

Let us denote the R matrix R(u — v) (||) considered as an operator V(j2u,) ® Y(H) — ► V(H) ® V(Uu) by 
R(Du,Dv)- The Yang-Baxter equation reads 

(1 <8 R(Ua,Uv))(R{aa,Uw) <8 1)(1 <8 ^(H.Ou,)) = (£(□»,□«„) <8 1)(1 <8 <8 1) 

: F(L\) (8 V(n„) (8 Vftiu) V^,) (8 V^) ® V(Du). 

We put (Fig. 0) 

i?(n ui ® □„» »•••»□«»,□«) := (^(a 11 ,a)) 1 ' 2 ( J R(Du 2 ,a)) 2,3 ---(«(o„ fe ,a)) fe ' fe+1 (n) 

: F(n ui <8D„ 2 (8 J -> V(D v )®V(n ul ®D„ 2 <8 ••• <8LT„J. 



Fig. 10 



The superscripts show the tensor components on which each operator acts. The resulting operator (11) is the 
operator which sends the rightmost tensor component to the left. We also put (Fig. O) 



RiUut ® • • • <8 DufcjUi <8 • • • <8 0u,) 

:= (^(Q,, ® •••»□„, ,rx ( )y- fc+ '- 1;fc+ ' ••• 

(£(□*„ <8 • ■ ■®0 Uk ,0 V2 )) 2 - k+1 -' k+2 (R(Bu 1 ® • ■ ■ <8 □ Ufc ,[] Dl )) 1 "' fci * +1 
: y^ ® • • • (giDuJ ® VfH,! (8 • • • ®DuJ -> V(Oui ® ■ ■ • <8QuJ <8 



(12) 



>□«,)• 



[ Fig. 11 ] 

For = 1, • • • , n, let l k be the Young diagram of vertical k boxes Q ( l 1 = )■ Then the fusion operator by 
Cherednik associated with l k is given by the comodule map corresponding to "the half twist" (Fig. |l2|) 

7r lfc := (R(a Ul ,D U2 )) k ' 1 ' k ■■■ 

(R(D U1 <8 <8 • • • <8 D Ufc „ 2 ,U Uk _ 1 )) 2 '" k ~ uk (R(U Ul <8 D„ 2 <8 • • • <8 □« fe _ 1 ,rJ u J) 1 "' fe_1;fe 
: y(n u -(fc-i)ft <8 •••<8D„) -» V(tlu®-"®O u _(fc_i) R ) (13) 

where the spectral parameters are specialized as 

(ui, ■ ■ ■ , Uk) = (u - (k - l)h, ■ ■ ■ , u - h, u) (14) 



[ Fig. 12 ] 

so that the rank of the operator ir^ degenerates. Since the operator R(D u ,Dv) actually depends only in the 
difference u — v, the operator ttih does not depend on u. The reason why this specialization ( p^| ) of parameters 
gives rise to the degeneration is as follows. First of all, from (|^) we have 

R{-n^ jf = -A(-h)%, (15) 



showing lm(R(— h)) C A 2 (C"). Actually for a generic value of Ti one can show lm(R(—Ti)) = A 2 (C"). Similarly, 
in tt^ one finds a degenerated factor R(—h) for each adjacent pair of tensor factors. With some braid manip- 
ulation, this implies the antisymmetricity of the matrix element of Ttik with respect to the permutation of the 
outgoing indices. This shows Im7r^fe C A fc (C") and in fact Imir^ ~ A fc (C") for generic h. We put 

V(l k u ) := ir lk {V(j\_( k _ 1)h <g> • • • <g> □„)) C V(D U O • • • <g> Ou-(k-i)h)- 

The fused R matrix as the comodule map V(l k ) ® V(l k ) — > V^l* ) ® is now defined as the restriction 

of the composed R matrix fll^), 

R{\ k u , l*f) := £(□„ • • • (g)D u _ (fc _i )ft ,a, ® • • ■ ® □ D -(k'-i)»)|y(i*)®vxi*') • 

We also define the fused L operator L(l^)|, € A(R) to be the matrix element of the comodule structure map. 
That is, choosing a basis {e 7 }/ C V(l*), we write 

l£ = ^(1*) 9 e J -> ^ e '' ® e V(l*) ® 

7' 

or equivalently, if we denote the dual basis to {e 1 } C V'(l^) by {e/} C V(lJ)*, we put 

:=<e I ,,l h u (e I )>. 

Since ~ A^F as a vector space, we can choose the basis elements in V(1JJ) as 

e 1 := sgn(a)e l ^ r > ® • • • <g> e v < fc > (J = • • • , i k \ix < ■ ■ ■ < i k } C {1, • • • , n}), 

<J6S(fc) 

where S(k) denotes the symmetric group of order k. Let {e^} C V* be the dual basis for {e 4 } C V, then for the 
above basis {e 7 } C V(l k ) the dual basis is just given by 

e/ := e 4l <g> ••• ® eij y(1 k) (7 = {ii, • • • < • • • < i k } C {1, • • -,n}). 

With respect to this choice of basis, the fused L operator is given by 

L(l k u )i = <e I ,,l k u (e I )> 

= < eg ® ■■■8e i ;,l^( ^ sgn(cr)e l "< 1 > <g> • • • <g> e l "< fc > ) > 

<JGS(fc) 

= < e i; ®---®e^,D u O---®n„_ (fe _ 1)fi ( ^ sgn(CT)e i "< 1 ) <g> • • • ® e M*> ) > 

crGS(fe) 

= J]] sgn(cr) < e^LUe^) > ••■ < , Q tl _ (fc _ 1)ft (e l -( fc ) ) > 

= J] S gn(a)L(u)^---L(u-(k-l)h)l w (16) 

trGS(fc) 

where 7 = {i\ < ■ ■ ■ < ik} and I' — < ■ • • < z' fc }. Applying the representation (|^) to this element we obtain 
a formula for the fused L operator acting on the space of functions on h*, 

A{R) 3 ml) 1 ,, ~ L(c\l k u Y r := Yl ^{a)L{c\u)\^ . . . L( c | u - (fc - l)*)^*' e End(0(h*)). (17) 

<JGS(fc) 



From the construction, they enjoy the relation ( |10[) . 
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3.3 Fused L operators from fused intertwining vectors 

The basic L operator (^) is defined as the composition of two quantities, the intertwining vectors. The fused L 
operator ( |l7| ) also have this factorized nature as we explain below. 

First we introduce "the space of paths" and their fusions to formulate the face weights as linear maps on 
appropriate vector spaces. We put 

•p/i-AP .~ / C : M = A + Tiei for some i, 
^ A 1 : otherwise. 

(Fig.|l3|) and denote by the basis of the one dimensional space V(U)x when n = A + hei for some i, and 
otherwise we set = 0. 

[ Fig. 13 ] 

For each u e C we consider the copy V(nu)\ of V(n)x and define (Fig. |lj) 

V (D U1 <g> ■ ■ ■ <g> Du J : = V (D M1 <8> • • • ® Uu J a ■ 



[ Fig. 14 ] 



These are the space of "admissible paths" in [ JKMO], in which the face weight and their fusions are formulated 
as linear operators on these spaces as follows. For g V(Ou)\ and e V(Ou)^ we put 



A u — v v 



thereby define the face operator (Fig|T|) VT(Du,H) : ^(Du ®0«) -> P(H <»□«), P(Du OD^a -> ^(H ®D«)a- 

[ Fig. 15 ] 

With these definitions, the Yang-Baxter equation of face type reads as follows (Fig.|l6|): 

(i ® ^{a u ,n v ))(w(n u ,n w ) ® i)(i ® w^n^n™)) = (vf(d„,d„,) ® i)(i <g> i4> (□„, d™))0^ (□„,□„) ® 1) 

^(□uOHODuOa -> T'COu, ® □«, ® CDa- ( 18 ) 



[ Fig. 16 ] 

The fusion procedure works for the face weight as well. We put 

W(D U1 ®Du 2 ® •••®Ott fc ,D ) := (^(□„ 1 ,^)) 1 ^(V^(D U2 ,D t ,))^ 3 ---(W r ([L fc ,D t ,))' £ - fc+1 



V{D U1 ®Du 2 



> □„„ (»□«)-> pen ® □«! Du 2 



(19) 



^(□u! ® • • ■ ® ® • • • ® n, ) 

:= (^(du, ® • • • ® □n fc ,D l , i )) fe - fe+ ^ 1;fc+i • • • 

W(D U1 ® • • • ® Uu k ,U V2 )f- k+1 * +2 {W{Ua 1 ® • 
: POXi ® • • • <S> ® □ Dl <S> • • • ® Uv J ->■ ^([l,! 



'□« fc ,D-ui)) 1 '" fe;fe+1 

• • ® Ouj ® [In ® • • 



(20) 



»Oui). 



and the fusion operator for l k is given by (Fig.|l7j) 
Hi* := {W(n Ul ,Bu 2 )) k ~ 1 ' k ■ ■ ■ 

(wiOa, ® n U2 <g> • • • <g> □„,, _ 2 , o Uk _ 1 )) i '" k ~ 1,k {w(u ul ® d„ 2 ® • • • ® , n Uk )) 1 "' k ~ lik 

■ P(Pu-(fc-l)R »•••»□„)-» P(du ® • • • ® □»-(*-!)») (21) 

with the same specialization (H) of (it], • • ■ , 

[ Fig. 17 ] 

Like 71]*, il]fc does not depend on it. We denote the image of this operator by V(l k ) := {V(n u -(k-i)h <S> 
• • • <8> Du)) C P(0« ® • • ■ ® □u-(fc-i)ft)- For J = {ji < • • • < j k } put ej := ej x H h ej d . We have 

p(l*)= p(l*)*+* e ', whereP(l fc )^ :^(D»®-®D„-( Mi )r'nP(l fc ). (22) 

A, JC{1, ■■■,«}, \J\=k 



The property similar to fllg) for the face weight with specialized spectral parameter — h implies that each 



subspace V{l k )\ +hKj (J C {1, ■ • • , n}, | J| = k) is one dimensional and spanned by the antisymmetric tensor 



e(1 fc)A+^ := sg n(a)e A (j ff(1) , • • • , j a{k) ) G V(l k ) x x +M ' (23) 

aeS(k) 

where 

If we denote the dual basis for e* +? ^ e V(Bu) X+M * by e *^ +fi ^ e (p(D«) A +agi )*, the dual basis for e(l k ) x x +Mj 
is given by 

e U JA — e AUl> ; JfeJ • e A ™ A+Re-jj 59 59 e A+fte^ +~+fte Sfc _ 1 h(l k )^ + J ' ' 

To define fused intertwining vectors, we formulate the intertwining vectors as linear maps as follows. 
Outgoing intertwining vectors (FigJT^): 



<t>(Uu) ■ V(Uu) -> V(p„), := ^ e J • 0(LV 



X3- 



[ Fig. 18 ] 

Incoming intertwining vectors (Fig.n9): 



[ Fig. 19 ] 

The intertwining property for these vectors read respectively as follows. 

%,a)<^(ri)®^) = ^>{Uu) ^ <P(Uu)w(Bu,Uv) ■■ p(Pu®o«) -> v(n, ®o„), (25) 
<^(ii)(»0(n u )^(n U! ^) = w r ([] u ,nv)0(Q 1 )®0(ri) : F(n u orx) ->p(a,8i],). (26) 

Putting 

</>(□«! ® • • • OOuJ := ^(Oux) <S> • • • ®0(DnJ, ^(dui ® • • • ®0«J := ^(DuJ ® •■ • O </>(DuJ, 
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[ Fig. 20 ] 



the relations (^5|), (|2(j) respectively imply 

7r lfc 0(Q Ul <g> • • • ®OuJ = <t>(Bu k ® • • ■ ODuJIl!*,, 0(Q Ufc ® • • ■ OO^Hi* = n lfe 0(D Ul <g> • • ■ ®DuJ (27) 

(Figj20|), where (ui, • • • , tifc) is as in (|l4|). Now we can define the fused intertwining vectors (outgoing and 
incoming) as the maps 

0(1*) := 0(tL fc ® ■ • • ® □ux)It'(iS) : ^(15) — ^(lj), 
:= 0(Du» ® • • • ® Du, ) | v(i») = - 

We also put 

0(l5)S == 0(l5)|p(i*)S;. ^(!«)a :=proi^(l*) 

where proj^ : V(l k ) — > P(l fc );( is the projection onto the direct summand in (p2|). ( proj^ — unless 
/i = A + he j for some J. ) 

The factorized nature for our fused L operator (|l^) is now stated as follows, 

Proposition 1 For J = {ji < ■ ■ ■ < jk}, let ej := ej x + ■ ■ ■ + £j d Tjf(\) := /(A + Tie,/). We have 

(L(c\l k u )j,f) (A) - E 0(l«+c R )A +tej /' ° 0(l«)A +fi£ " J ' I (^/)(A). (28) 

JC{l,---,n},|J|=fe 

//ere o stands for the composition of : V(l„) — > P(l^) and 0(1«+ c r) : ^i^u+ch) ~ * ^O-u) v ^ a ^ e obvious 

identification V(l^) — 'P{^u+ch) > super/subscripts 1,1' denote their matrix elements. 

Proof of this proposition is rather straightforward from the definition and omitted. □ 

4 The commuting family 

Now we are in a position to ask the question in Section |l|: "what kind of operators arise as the traces of these 
L(c\Y u ) 's?". We shall consider the case Y = l d , the vertical d boxes case. Then L(e|l„) is a matrix of size 
dim A d C n whose matrix elements are difference operators. 

Theorem 2 (Main Theorem) Let M d (c\u) := Tr ld £(c|l^), d = 1, • • • ,n. We have 

A,e s -e t > +- 
. K < A, l s - e t >) 

JC{l,-,n},|I|=<i \s<£l,tel y ' 

where 0(u) is the Jacobi theta function (Qj and Tj stands for the h-shift operator: 

(Tf/)(A):=/(A + ^), :=[]??■ 

iei 

These operators form a commuting difference system: 

[M d {c\u),M d ,(c\v)]=0. 

Remark. It is an important point that these operators obviously commute, as we mentioned in the intro- 
duction. This is because the extended "RLL=LLR" relation ( |l0| ) can be rewritten as 

R(Y U , Y^)L(c\Y u ) ® L{c\Yl)R{Y v XT X = L(<M) ® L(c\Y u ), 

and then taking the trace simply gives 

M(c\Y u )M(c\y:) = M{c\Yl)M{c\Y u ) 

where M(c\Y u ) :— Try ry) L(c\Y u ) . This simple argument and the resulting operators, "the commuting transfer 



matriecs", was effectively used in Baxter's analysis of the spin chain models |Bax71|. Thus our result can be 
said as an ideology: 
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commuting transfer matrices = commuting difference system. 



Proof of Theorem [2|. From Proposition [l], the operator is of the form 
M„(c|«)/(A) = L{c\lX +h ^/{T r }f){\) 

I, JC{\,— ,71}, \I\=\J\=d 
JC{1,— ,n},\J\=d 

By the cyclicity property of the trace, the coefficient for Tj can be written as 
Tr v(ld) (Hl d u+ch ) X x +h *-> o 0(1^+^ : - p(l<*)^ - 

Now taking the trace is much easier because r p(l d }* +hej i§ one dimensional. We take its basis (^3|) and the dual 
(pif). Moreover, we note 

^< +c n)l +Mj {ni d )l +Mj ) = m+cn ® • • • ® □ u+cS -( d - l) ,)^ ft ^(p(i d )^ fiEj ) 

C FCDu+cR <g> • • • <g> Oa+ch-(d-l)h) V® d . 

The last isomorphism as vector space is induced by the obvious identification V"(Du) ~ V. Then we have (Fig.|2l|) 



TV. ^ (KO\ +hlJ ° <f>(l d u + cH)x +m ■ V(l d ) x+hSj -> V(l d ) -> V{l d ) x x +M ') . (29) 



Tr v(ld) x+ 



7 ® • • • ® □„-(d-l)ft)A +R "' ? ° 0(tiu+di ® • • • ® []„+cR-Cd-l)R)A 4 



v{i d ) x . +hi 



'A ' ' ' V 1 )\ 

= < e*{l k ) x + n ^'^{Uu ® • • ■ ^-(d-Dfc)^ o </>([]„ +c fi ® • • • ^Q u +cn-(d-i)h)l +m (e(l k ) X+h ' ej ) > 
= E s S n ( a ) < e *\(h,---Jd), 

oGS(d) 

4>{Ua ® ■ ■ ■®B u -( d -i)h)x +hSj ° HOu+cn ® • • • ®Uu+ch-(d-i)h)l +hEj (e X {j<y { l), ■ ■ -Jo-id))) > 

creS(d) 

<x6S(d) fc=l \m=l / 

[ Fig. 21 ] 

We are in a position to use the following formula 



^)£ + ^'">(u+c?J)^ +, '' 

rn—l 



n 



U i) 



Lh 



■ 1 1 1 



(3 + 1) 



th 



(S- < AMi-i >) fcC 2 **- < A,e" >) < M> e i+X >) 



< Mi ei-i >) On(^~ < A, >) 6 n (%- < M , e i+1 >) 



•!**(£-< am* >)L, fc , =1 ,., n 

f + u+ < fi, ej>-< A, e, >) 



n 



f + < /i,ey > - < A,Ci >) 



(< fl,6j> > - < fl,6j >) 

which follows from the definition of 4> (||) and the determinant formula of Vandermonde type 

6j(u k ) 



det 



= ( 1)n _ 1 g(Sj"j) T-r 0{uk - Uj) 



For convenience, for any statement P let 

Y P := 1 if P is true, F P := if P is false. 

We have 



(31) 



(32) 



A+fte^, H h^Cj 



, A+ftej H hTiej 



6{f+u-(k- l)h+ < A + he h + ■ ■ ■ + he jk _ t , e Jk > - < A + ^ + ■ ■ ■ + Tie 



(tt-(fc- i)n) 



'6+< A + fte,- 



> - < A + ^(i) 



fie 



J<r(fc-1) I C J<r(fc) 



iVifcSj'<n ^ A + + ' ' ' + ae ^~i , £j' > - < A + fte^ + • • • + fie Jfe _ 1 , ej fc >) 
0(£ + u - (A - l)h+ < A, e- - e Uh) >) ^ 0(f + < A, e> - e jCT(fe) > +»> e{jl ,.. Jfc _ l} ) 



0(u- (fc- 
0(f +u-(k-l)h+<\,e jk -e Uk) >) 



n 



6>(< \,e f -e jk > +KY j , &{ju ... tjk _ l} ) 



n 

r^fe,l<r<d 



6{u-{k-l)h) 

(f + < - CjVw > +^r<fe) t-t 0(f + < A,e> - e Uk) >) 



< A, e 7r — e,-. > +7ii^ 



<fej 



n 



0(< A, £j v - e jk >) 



Therefore 



■ ch 



+ U-(k- 1)7?,+ < A, e 



Jk e jcr(k) 



k=l 



9(u-(k-l)h) 

;f +< A,e ir -e jV(W > +fir r<fc ) 



n ' 

rjtk,l<r<d 



(< A, e jV - e Jfc > +?iY r<fc ) 



n 



' ch 



+ < A,e 



J C J<r(fc) 



>) 



e(< a, ty - ej k >) 



Eaes(d) sgn(^) nLi Ilr=i (f + < A ' e> - ^ w > + «k,r(« - (* - m) 



7£S(d) 



l<r<k<d ' 



(< A, e jr - e jk > +h) II 



l<fc<r<d ' 



(< X,e jr - e 3k >) 



n 



0(^+<A,e>- £ - 3 >) 
0(< A, e r - ei >) 



jeJ,f<£J 

det [nli #(f + < A, e jr - e jy > +KY r<k + 5 k>r (u - (k - l)h)) 



yrd-1 
llr=0 



fi)n Kr<k<d 0(< A.e^-ej, >+fi)I[ 



l<fc<r<d ' 



(< A, e ir - e jk >) 
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n 



0(%+<\,e r -e j >) 
(<-Mj'-^>) 



(33) 



To end the proof, we are led to show the following Lemma |l| below, which is quite interesting itself. According 
to the lemma, we see the unwanted factor in (|33]) cancels as desired. QED. □ 



Lemma 1 Recall Y r<s := 1 if r < s holds, Y r<s := otherwise. The following formula holds: 

d 



det 



Y[ {n r - \ s > + KY r<s + S rtS (u - {s - l)ft)) 



(34) 



s,s / =l,"-,d 



d-1 



l<s<s'<d 



Here the variables {A s },{/i s } are arbitrary, i.e. neither the conditions y] , A a — nor ^ s=1 ^ s — are 
necessary. Proof of this formula uses the induction on d and will be given in the Appendix. The h = case of 
([mJ) is easily transformed into the Cauchy type determinant formula 



dct 



9(fi s - \ s > + u) 
6(n s - \ s ,)0(u) 



8(u + Y,LM - V)) Ui<s<s'< d 8(»s - /v)0(A* - A s ) 



s,s' — 1, ■ ■ -,d 



6{u) 



n s ,y=i,., d % s -A,o 



(35) 



which is also known as the genus 1 case of Fay's trisecant formula [Fay]. But I do not know whether (|3J) is 
previously known or not. 

It is also interesting to remark that this h = case was quite relevant in Q although his approach for the 
commuting system is different from ours. 



5 Relations to other approaches for the system 
5.1 Ruijsenaars' operators 



Proposition 2 Put g 



Define a function <£> on h* by 



oo oo 



*(a) : = n d + ( Zk /z k ,), d+{z) : = n n T 



1 - zq m+1 p k 1 - Z - l q m - g p k+l 
— zq mJrgJrl p k 1 — z~ 1 q nl p k + 1 



(36) 



where p = Cjxp27riT, q = exp 2-7ri?i(|g| < 1) and Zj := cxp 2ni < A, e,j >. Then the conjugation by the square root 



'ft 

yields u Ruijsenaars ' [[rJ commuting operat 



9{v) 



®- 1/2 M d (c\u)$ 1/2 



ic{ 



V ( TT l e &+<^ts-e~ \ ( n l 9(fh+<X,e t -e~ \ 
^ 11 V ^(<A,e s - ei >) MM H V d(<\et-e s >) 

l,-,n},|/|=c2 \s<£LteI » v ' / \s£l,tel V v ' / 



Proof, [cf. ||Die|] The function d + (|36|) satisfies 



1 - zqp k 1 — (zq) 1 p k + 1 



d+(qz) iJL i _ zq i+gpk 1 _ ( zq i+ g yi p k+l 



as long as it makes sense 
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and consequently $ satisfies 

$(z) $(zi,---,z„ 



j (1 - ^)(i - |g °° a - (g ) P fc )(i ^ (grvxi - (g - (f g 9 )- V) ' 

11 j (1 _ SffiL g ,)(i _ |) 11 (1 _ (£^)p*)(l - (^q 9 )"V)(l - apfc)(i - (|)-V) 
6(H + \ij)6{gh + 



where Xij —< A, — e"j >. From this property, for general / C {1, ■ • ■ , n} we have 

_ -q e(n + Xi tj )e(gn + x jti ) 



T?${z) ~ Jj^ 9(gh + h + .\,.,;fl:A ; ... ■ 



Therefore, 



( U + M)\ 1 



(V^M^u)^ 2 /) (A) 



£ f II eiX 8(\ + ) h) ) ■ (*" /2 ( T ^) 1/2 ) W(2?/)(A) 



^(A,,i + 









0(A^ + 




^(Aj,-. 





0(A A < + 




0(X 3 ,, 





e{Xj,i + 





1/2 

e 1 n 

|/|=d \iei,j^ 

\ 1/2 / \ 1/2 



1/2 / ^ 1/2 

0{gh + X tJ ) 



□ 

Remark. Put := d + (z/q),S := Yik^k' d(zk/zk>) and define /(A) := /(—A). Consider the inner product 
defined by (/, g)' :— [fgS]i, where [ ]i stands for the constant term in A. Then generalizing the trigonometric 
(Macdonald's) case, one can directly verify that the operators {Md(c\u)} are formally self-adjoint with respect 
to this inner product. 



5.2 Krichever's Lax matrix 

Moser system (for the proof see (|48|) in Subsection I 



Put Md := e { u +£d h -j Md. In the limit % — > 0, one can check that our system degenerates to the elliptic Calogero- 



(-2M 2 + M\ - 2Mi + nj ^ H, (37) 

H = AC/n ° ( E £ - K + 1) Eoos'H^) ] ° A ~ c/n - 

i— 1 J i<j 
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Here A(A) := Yl%<j Note that the Jacobi 9 function and the Weierstrass sigma function with quasi-period 

l,r are in the relation a(u) = e hu 9(u)/0'(O), h = —2n^—lJpp log 77(7-) so that we have p(u) — — (log<r(u))" = 
— (log f?(u))" + h [CH|. It is well known that this system is integrable and admits a Lax matrix formalism jKiJ]. 
Here we shall derive the Lax matrix from our L-operator. 

It turns out that the L-operator (||) should be slightly modified for our purpose. We conjugate it by the 
matrix formed by the intertwining vectors: 

n 

L(c\u))^L(c\u% := Hu) X x +M ^Hu) X x +M ^L(c\uY r 
Using the property (||), this conjugated matrix takes the simple form (Fig. p2| ) 

n 

L{c\u)) = £ ^u)l + ^ k d>(u + ch)\ +m * k Tt 
fc=i 

Moreover, for the summation in the right hand side we can apply the formula (|3l|): 

[ Fig. 22 ] 



f +« + A i)f ) yr9(f +A M ) 



fe=l 



so that we have the formula 



0(A fej ) 



L(c|«)5 = 



n 



9(Xk, 3 ) 



rph 



This is our version of the Lax matrix suggested in |r| . From this it is easy to see: 



H—>0 



L(c|u)J^a}. 

To recover Krichever's Lax matrix, we look at the coefficient for h 1 in L. We have 



d ( 9{% +u + X j , i )_ ~6(% + A M ) 



*(«) 



d (9 (f + u + X j , i ) n g(f +A fc|i ) 



an 



J at. 



T ■ 



h=0 



c 
n 



+ + g (" + A^) 



c /0'(u) 5 



^ ^4 + (1 _ ^'(o) nfe ^ e(Afc < 



0(A. 



rt \9(u) de, 



■ lOg A A + _ I + (1 _ gj) w\ N TT fl/x 



(38) 



(39) 



If we further modify this matrix by 1) the conjugation by A(A) c / n and by 2) the similarity transformation by 
the diagonal matrix diag ( Jlfc^i 0(Afc,») ) , we obtain 



Proposition 3 



A c/l 



?i=0 



c (9» 5 



n 0(u) 3ej 



+ ^ + (1 - 



c g(it + A&QgCO) 
n e(u)e(Aj,i) 



AT(c|u) 
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One sees that this matrix K(c\u) essentially coincides with Krichever's Lax matrix modulo the scalar matrix 
addition £ T$i<L 

In the classical case, the so-called "dynamical r-structure" for Krichever's Lax matrix is obtained by Sklyanin 



[393]. From our viewpoint, this can be easily extended to the quantum case and the dynamical r-matrix can 
be essentially understood as the first derivative of the face weight at U = 0. Thus the dynamical variable 
dependence of the dynamical r-matrix quite naturally arise since the face weight depend on the variable A 6 h*. 
We will treat these points in a seperate paper. 

It should be also mentioned that the m atrix L(c\u) gives an representation of the quadratic algebra in pq| 
after a suitable formulation. It is known [BD] that the elliptic extension of the vertex type solution to the 
Yang-Baxter equation only exist for the type A case (i.e. the Belavi n R-matrix), while the face type solutions 
(face weights) are known to exi st fo r the other types as well JMQ2| . This may suggest that, to generalize the 
Olshanetsky-Perelomov system] OP] with the symmetry other than type A, searching the analogue of the L(c\u) 
matrix might be useful (cf. [ABB]). 



5.3 The generating function and the differential limit 

Conceptually, our commuting operators {M d (c\u)}2 = i are generated by the L-operator L(c\u) — [£(c|u)*-] by 
fusing and taking the trace. In this sence L(c\u) is the generating function of the system. More explicitly, we 
can state the relation in the following familiar form: 



Theorem 3 



J2{-t) n - d M d (c\u) = :det[L(c|u) - t] 



i=0 



where we put M (c\u) := id and the normal product :: is defined by putting the difference /differential operators 
to the right: e.g. :f(\)Tfg(\)T r }: := f{\)g(\)T r fT 1 } for any function f,g in A. 



Proof. As in the trigonometric (Macdonald's) case, we can rewrite our operator as follows: 

G(u) 

I\=d 



M d (c\u) = 



6{u) A 

where A(A) = Y\ i< j ^(A? — I* 1 view of the determinant formula (|32|), we have 
8(u + T-" h A n . _ T-* h det[9j(% - X. t 



Therefore, 



6{u) 



J2(-t) n - d M d (c\u) 

d=Q 



■(A) 



detie^-Xi)} 



for 7C {!,••-, n},\I\ = d. 



EH)"- d E 
\i\=d 



'Tj -"det^-Ai 



d=0 



det^-Ai)] 



TV 



EH) 



n _, f l =det[(7;- 



' hYieI 0^^. — \i))Tf YieI ] 



^ctfc^-A^-t-fl^-A.)]: 
det^-A*)] 

:det [9 3 (^ - A,)Tf — t • 6j(~ — A,)] ■ det ~ A,)] _1 : 



:det 

:det[L(c|w) -t] 



(40) 
(41) 
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Remark. Let 1 := J2i=i £ i- Then [Md(c\u),T^ ] — for any a. This ensures that the replacement of 
rpah _ ex p fl ^^_ ex p a ^^_ m f ormu l a f Md{c\u) will not affect the commuting nature of the family. 
Thus we can and will abuse the symbol T°" h (resp. di) for these two meanings, expahei and expaTie^ (resp. JjL- 
and £§-), in what follows. □ 

The numerator in ( [40| ) can be regarded as the elliptic difference analogue of Jiro Sekiguchi's generating 
function for the trigonometric differential case pej ]. 

In (pl|), there is another option, i.e. multipling det [0j(— — A;)] 1 from the left. As is easily seen, the result 

is 



J2(-t) n - d M d (c\u) = :det 



d=0 



:det[£(c|«) -*]:, (42) 



,fc=i, 



where L is given by (|3q). Using this generating function, we can derive the commuting operators in the 
differential limit Ti — + 0. Since L id (|39|), ^(L — 1) n —> Q -^Lfi{c\u) =: L' (c\u), we have 



:det[L h (c\u) - (1 + fit) ]: = :det[n(L&(c|u) - t) + o(fi)]: R ~° fi n :det[i' (c|«) - t]: + o(fi"). 
This leads us to define the differential operators {Dd(c\u)} by 

L h {c\u) - (1 + U) 



lira :dct 



:det[L(,(c|u) - t]: =: ^(-t) n - d ^(c| W ). 



(43) 



d=0 



In terms of Krichever's Lax matrix (Prop. 3), we have an equivalent definition 

n 

A c/n o :det[K(c\u) - t]: o A~ c/ " = ^(-i)™- d D d (c|u). 



These formula actually define a commuting family [Dd(c\u), Dd'(c\v)] = 0, because the left hand sides are just 
the summation of the commuting operators {Mj}- Obviously, Dq = id. The higher operators are given by: 



:det[L' Q -t}: = ^H)""" ]T : det[Z' (c|<] w : 

IC{l,—,n},\I\=d 



d=0 



J2(-tr- d J2 E detK( W ,A)^, /eJ .^, 



2=0 



|J|=d JC/ 



where := Jlie/ ^ an< ^ 

C' {u,X)) := Ac fi (u,A,A)^ 

Since CS(«,A)5 = gf-C (u, A, ^) 



C n (u, A,/x)5 := -A» S7 -f^ 1 ± JT 7 (of.©). 



X — /.L 



we can apply the formula (B5h as follows: 



detK(u)^, /eJ = (-f)l J ldet 



gf-Co(w, A, /i) 



X—fi 



(~0 Jl (ni eJ ^)det[C ( W) A, M )] i , eJ 



X—fi 



(-i ) |J| (n 



e(«)A J (p)A J (- At ) 



(44) 



(--) 



v ( -i)i^^M^(A). 



X—fi 



J=J'UJ 



Here A, 7 (A) := Ujj'eJ,j<j' 9 ( X i ~ V) ■ Summarizing, (note that = ^A for J" C J C {1, • • • ,n}) 
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Proposition 4 Put (—^d) J := Iljej — ^ e differential limit fi — > 0, i/ie commuting operators {Dd(c\u)}^ = 
are given by the formula 



D d (c\u) 



E <-,,M-.!7,.,^«A,.(-n*)™'- 

|/|=ci J'UJ"GI 



I\J'=:I' 



d J A 



Ehi'Vw E E 



fc=0 



J'|=d-fc J"C/' 



■2fl 

c 



IV" 



(45) 



Since the functions {^-{u)}k=o,---, n are linearly independent, we can consider the coefficient for (— 

it is for d < k, and 

E E^W-(-^)' V =:^-*] (46) 



in (— f )™ -Dd(u) above. From (^5|), it is for d < k, and 

d J A . / // n ' ' 

|7|=d-fc Jcl 

for d > fe; it depends only on the difference d — k. For example, 



up] = 2(--ft)(--a i ) + -L_(_-a J ) + ^_(--a f )+ ■ 



«<3 



A 



These are the operators appeared in [ Dc ] in the trigonometric case: 

Corollary to Proposition ^] For d= 1, ■ • • , n, the operators D[d] commute with each other. □ 

To end this section, let us prove the previous formula (]37|). As is easily seen, the "hamiltonian" H and these 
operators are related by H = D[l] 2 — 2D[2]. On the other hand, 



Proposition 5 



D 



fe'=0 



n - k'\ {-^) k d d n - k M k ,(c) 
n-d) k\ (d-k)\ 



Proof. By the defintion of D[d — k], we have 



d>k 



From (|43|) and Theorem ||, this is equal to 



fc'=C 



7i n 0(u) 

n-fc' 



fe'=0 



d=0 



d 



M*(c) 



dt-d—n 



oo /ck'%\k 



t a Ti 



E 

\fc=0 



fc! 



(u) M k ,{c) 



(~) EMr fc 'E n 7> d E 

fc'=0 d=0 V 7 fc=0 



fc! 9 w (rc-d-fc)! 



n=o 



Rearranging the last summation to obtain the coefficient for (— t) n ~ - — 4 («), we get the assertion. □ 
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Some particular cases of this proposition are: 
D[l]=D[l-0] = 
D[2] = D[2 - 0] = 



( n o 2 )M2 - (V) Af " _ Mg - (n- 1)M{' 



where ' = dn\n=o- Now the formula ([37]) can be verified as follows 
Ml - 2Mi - 2M 2 + n 



A - n ^ 2 , 9 fo~ - n) - (Af 2 - g)j 

^ M? + 2 {n - 1) f-^&D [ lf-2D [ 2] = H. 



(47) 



(48) 



In the last line we have used the l'Hospital rule (since (n — 1)M[ — M' 2 = 0, we are to take the second derivative). 



6 The invariant subspace spanned by theta functions 

In this section, we would like to investigate a certain finite-dimensional invariant subspace for Md(c\u) when c 
is a nonnegative integer 

Let us consider the space of level I theta functions on the weight space h* , 



Th, 



f : h* C 



f(X + a) = .f(X), 

/(A + ar) = /(A) exp[-27riZ((A, a) + (a, a)r/2)} 



(Va S Q) 



(49) 



where Q = Q(A n _x) stands for the root lattice, and also let Thf^ be the subspace consisting of symmetric 

group action invariants. This space is of dimension and actually spanned by the level I A^lj -characters. 

For j = 0, • • • , n — 1 let Aj be the classical part of the j-th fundamental weight for the type A^_ x root system 
[ Kac and let 

XjW ■= exp27r\/^T 



Then more precisely we have 



Th 



In [H2|, we have shown that Thf is an A(R) submodule in 0(h 



S(n) 



e 

0<3i<-<j'j<n-l 



{\,V) 

CXn 



(u, v)t 



■Xji 



(50) 



Theorem 4 ([H2|) When the parameter c is set to a nonnegative integer I, we have 

L(Z|u)}T/if (n) c Thf (n) 
for any u € C and i, j = !,••• ,n and therefore 



M d (l\u)Th 



S(n) 



C Th 



S(n) 



(To be precise, the operators L(/|u)j as functions in u have poles at u E Z + rZ and therefore the space 
Thf^ should be extended by tensoring 0(C), the field of meromorphic functions in u S C.) □ 

Here we wish to establish a structure theorem for this subspace as an A(R)- module. 

To describe the result, let us give some definitions. First we introduce an A{R) module which is isomorphic 
to the A(R)- comodule V(U Vl <8> • • • ® H,) as a vector space. To avoid confusion let us denote the space by 
V(Dvi ® ••• ® Duj)- The A(R) module structure for this space is given by [ FRTfl 



L(u)j,e»* = E,v e>' ■ R{u - v)*A, for e^" e V(Du) 



19 



and in general we tensor this structure via the co-multiplication 

L{u)\,^ Yl L(u)} 1 ^---^L(u)^- 1 €A(R) & (51) 

so that we have 

L{u)i,v:=R{n u ,U 0l ®---®n Ul ) i i ,v for any v G V(D Vl ® • • • ® Q v , ). (52) 

Here JiflHajD,,! ® • • -SSH,)!/ stands for the operator on V(Dtn ® • • defined by R(n u ,U Vl ® ■ ■ •®Di;,)e*®v = 

e *' ® (-R(Dm,[V ® • •■ ^DoJi'V). Recall that lm(R(-hj) C A 2 (C n ) C V(0« ® H-ft) ©■ For more general 
lengthy case, we have 

l-i 



X := £ Im^(-»)* , * +1 C V(0, • • ■ ® []„_ (l _i )R ) 



fc,fc+i 



m— 1 

and the subspace X is an A(ii)-submodule. This can be easily seen from the relation 

(k-l)th k-th . v 

i?(Du, IX ® • • • ® Dv-(k-i)H ® □ v—kh Qv ■ ■ ■ 

(fe-i)th fc-t/i 

= J?(-n) fc ' fc+1 i?(0«, H ® • • • <8> Ck-fcfi ® D«-(k-i)fi ® • • • ® , 
which is of course a consequence of the YBE. Now we define the quotient j4(i?)-module 

S l v := V{U„®---®n u - ( i-i) h )/l. 

By the property (|l^) we have 

■ ■ ■ (g> {e l (g> e j - e J (g> e l ) (g> ■ ■ ■ = mod X, 
and hence is isomorphic to the usual space of homogeneous polynomials of degree I in e 1 , • • • , e n ; 

S l v = C e?' 1 • • • e Ji , e jl ■ ■ ■ e n := e jl ® • • • ® e Ji mod X. (53) 

0<ji<-<ji<n-l 

Comparing fl5C| ) and (p3|), they are obviously isomorphic as vector spaces by the correspondence 

7 : e J1 • • • e Ji i-> • • • x i; . (54) 

Theorem 5 Let us modify the L-operator ^) in the case c = I, a positive integer, by the scalar function 
multiplication 

L °m---(li 9 -%f 1 )-mu)) (55) 



and consider the A(R)-module structure on Thf defined by L(u)j i— ► L°(l\u)j. Then the correspondence (]^| 
gives the isomorphism of A(R) -modules 



7 : S l Q ► Thf (the parameter v for S l v is set to 0). 



Remark. In the case I 



1, we have shown in [EE], Theorem 2(1) that 

0(h) 



6{u) 



(56) 



implying the isomorphism. This immediately explains why the modification ( |55| ) is natural. From d56j), it follows 
directly that Xj(j — 1) ' ' ' i n ) are eigenfunctions for Mi(l|u) (and actually for the system {Afd(l|u)}2 =1 ) with 
the same eigenvalue jfil J27=i R( u )lo- 
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Proof of Theorem [||. Based on the / = 1 case (see the remark above), the general case can be verified as 
follows. 



7 (L(u)*,(e*-..e*)) (A) 

7 (R(Ou, Do • • • (e 31 ® ■ ' ' ® e Ji modX)) (A) 

7 I E (i?(u + 0)^ 1 e-'' 1 )(»(i?(M + ?i)^e j2 )(8)---®(i?(u+G-l)n)^- 1 e j! )modx] (A) 

,ii_i=l 



= 7 E (iJ(«+o) < i y i )(ii(«+R);^ a )... (Ji(u+(i-i)a){}-v«) Wa) 



n 



E + °)i 1 Xji)( i °( 1 l w + ft^Xjs) • • • + - l)ft)i' l-1 Xj'i) ) (A) (57) 

>»i-i 

it + s/i) 



1-1 



X E E (<K")A +r " fcl V(« + ^ +ri %M- + ^^(tl + 2^)^- 2 ) 

*ir">*l-i = l fei,-")fcj=l 

• • • {4>{u + l)hf x +Mk ' l '-^(u + in) X x +m 'i') ■ Xji(A + ?ie fel )Xi 2 (A + 7ie fc2 ) ■ ..^-.(A + fie*,) 
j-j- #(it + sTi) 

s=0 



9(h) 



x E ^Ma^ 1 ^^,^ ■••4 1 _ 1 ,fc>(M + "i)A +S£ ''i' ' Xji(A + ?ie fcl )---Xj ! (A + ?ie fc! ) (58) 
fci,---,fci=i 

LI TWT^ • E + ^)a + *% (Xi^ • • • Xi, ) (A + te*) 

s =o a[n) fe=l 

(L (I|u)J,(x* ' ■•*«)) ( A ) = (i°(^)W(e J1 • • • e Ji )) (A). 



[ Fig. 23 ] 

We used the I = 1 case in (|57j). The line ( f58| ) is the most essential step, where as shown in Figj2^ the situation 
allows us to employ the duality relation for the intertwining vectors □ 

Thus there is an interesting "representation theoretic" invariant subspaces for our operators. This space 
would be identified with the space of Weyl group invariant theta functions in [ EK2| , where they considered 



an affine analogue of Sutherland operator and its diagonalization. Note that, since our operators recover the 
elliptic Calogero-Moser system (p7|), our system can be regarded as the q- analogue of the critical level case of 
Etingof-Kirilov's theory. Thus the existence of n commuting operators Mk(c\u) is in a coincidence with the 
large center of the (completed) quantum affine enveloping algebra at the critical level. 



7 Discussion 

In this paper we realized Ruijsenaars' commuting difference system as commuting transfer matrices in the 
solvable lattice model point of view. In other word, this system can be considered as a function space realization 
(a Schrodinger picture) of the one-dimensional spin chain defined by the elliptic R-matrix as its Boltzmann 
weight. Regarded as a lattice model, this system is very trivial because it contains only one site of freedom. 
It is quite natural to think of the multi-sites case as well, that is, consider the (fused) L-operators and their 
traces acting on the tensor space T/if (n) <g> •• • <g> Th^ n) . On this function space we can consider the action of 



21 



the commuting transfer matrices again. The resulting system can be regarded as a difference version of the 
elliptic curve case of generalized Gaudin models of type A :jEl|, gY|, jFFRj ], pNfl , Q. In these papers Gaudin 
Hamiltonians are generalized upon a geometric setting, and the commuting Hamiltonians can be obtained by 
the pole expansion of the characteristic polynomial of the (quantized) Higgs field. Therefore the factorized 
L-operator or the difference Lax matrix L (Proposition 3) can be said as a difference analogue of the Higgs field. 
This suggests that the intertwining vectors, the building block for the factorized L-operator, would also allow 
some intrinsic definition (geometric as well as representation theoretic). The paper [KrZ| is in this direction 
with the geometric treatment of the intertwining vectors and consequently the factorized L-operator for Baxter's 
R-matrix. But still it seems that th e orig in of the weight space, or the dynamical parameter A, needs to be 
elucidated. Another related work is | C95|] , where Cherednik announced the construction of elliptic difference 
operators based on his Hecke algebra representation technique. Whether and how we can find the relationship 
with his thoeory and the present construction would be also an important problem. 
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A Proof of Lemma 1. 

Recall y r<s :=liffr<s and otherwise, and put 

d 

A W = A («0( u> Aj m)s s , ;= Y[ 6 ( Mr - A s , + TiY r<s + S r , s (u - (s - l)7i)) 



Then 



A{ d l 



s,s / =l,-",d 

that d — 1 case is true. 



is the matrix in problem. Lemma [l] ( |34| ) trivially holds for the case d = 1 ; suppose 



First we let /id = A s ' for some s' = !,■■■ ,d, and observe 



• A^ s , = 0(s < d), Aft = 9{u - (d - l)h) UrZi 0{»r - A s < + h) and 
. A<5 = Af- l) ■ 9(X S , - A s ) (1 < M < d - 1). 
Then the determinant is simply a product of A^p, and the remaining (d — 1) x (d — 1) determinant: 



,(<n 



rd-l 



d— 1 



det 



.4 



= dot 



s,t=l,— ,d 



A 



(d-1) 



s,t 



■ 6{u -{d- \)h) TT 0(/i r - a s , + h) TT 6{\ s , - a.) 

l,—,d-l;tjts' -M; J v- L 



r=l s.#s' 

By the induction hypothesis, the right hand side is equal to 

d d-l 

6{u + Y^&r ~ A r )) 0(u - sh) Y[ 9(\ s/ -\ s )6(h + n s -fi s >) 

r=l s=l l<s<s'<d 

as desired. That is, the formula in question holds for d special values /id — A s < (s' = 1, • • • , d) of fid- 

A standard complex analysis tells that, for an entire nonzero holomprphic function f(z) with the properties 
f(z + 1) = exp[-27r\/^LB]/(z), f(z + t) = exp[-27rV^T(C*z + C')]f{z) 

have C zeros z\ , • • • , zc with multiplicities inside a fundamental period parallelogram and 

C 



^ Zj =Bt+ C mod Z + Zt 
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holds. 

Consider the difference 

d-1 

- 9(u + y^(^ r - KY. 

,t=l,--,d 



det 



l<s<s'<d 

as a function in /id = z and denote it by f(z). Then it is easily checked that f(z + 1) = (— l) d f{z) and 

d 

f{z + r) = exp[27r\/=T(-u+(d-l)a)] J](-exp[27rx/=l(-Md + A 5 -r/2)])-/(^) 

/(*)■ 



exp 



d j 



Therefore, if / is not identically zero, we should have d zeros modulo the period lattice whose summation should 
be given by 



d d , A. ,r + l 



r + - - (- A s + d— h u - (d - l)fi) = A s - u + (d - l)h mod Z + Zr. 



On the other hand, the previous observation shows that we have d zeros Ai,---,A<j whose summation is 
J^^ =1 A s ^ Ss=i A s — u + (d — Since u and 7i are arbitrary, this is a contradiction. Thus / should 

be identically zero, proving the lemma. □ 
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Figure Captions 



R(u - v)*, = R(u - v)% 
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Fig.|l]. Matrix elements of the R-matrix. 
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The outgoing intertwining vector. 
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The incoming intertwining vector. 
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The duality relation. 
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The factorized L-operator 
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Fig.0. 

The intertwining property of the outgoing vectors. 
Dots represent summation indices. 
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Fig.0. 
Face fusion operator. 
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Calculation oTthe trace. 
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The matrix L(c|u) 
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